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ABSTRACT 

In this paper, we investigate the structure of condensation fronts from warm diffuse 
gas to cold neutral medium (CNM) under the plane parallel geometry. The solutions 
have two parameters, the pressure of the CNM and the mass flux across the transition 
front, and their ranges are much wider than previously thought. First, we consider 
the pressure range where the three phases, the CNM, the unstable phase, and the 
warm neutral medium, can coexist in the pressure equilibrium. In a wide range of the 
mass flux, we find solutions connecting the CNM and the unstable phase. Moreover, 
we find solutions in larger pressure range where there is only one thermal equilibrium 
state or the CNM. These solutions can be realized in shock-compressed regions that 
are promising sites of molecular cloud formation. We also find remarkable properties 
in our solutions. Heat conduction becomes less important with increasing mass flux, 
and the thickness of the transition layer is characterized by the cooling length instead 
of the Field length. 
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1 INTRODUCTION 

It is well known that in the interstellar medium (ISM), 
a clumpy low-temperature phase [cold neutral medium 
(CNM)] and a diffuse high-temperature phase [warm neu- 
tral medium (WNM)] can coexist in pressure equilib- 
rium as a result of the balance of radiative cooling and 
heati ng owing to external radiation fields and cosmic 
rays ((Field. Goldsmith fc Habind Il969l : IWolfire et all 1 19951 . 
120031 ) . The CNM of atomic gas is observed as HI cloud 
(n ~ 10 - 100 cm" 3 , T ~ 10 2 K). These two phases are 
thermally stable. On the other hand, the thermal instabil- 
ity (TI) arises in the temperature range between these two 
phases. Thus, the ISM in atomic phase can be interpreted 

as bistable fluid. 

Pioneering work on the TI has been done by [Field 

( 1965) who performed a linear analysis of the thermal equi- 
librium gas and derived a simple criter ion for the TI. Fo- 
cusing on a fluid element. iBalbusI (|l986l ) derived a criterion 
for the TI of thermal non-equilibrium gas. In the nonlin- 
ear evolution, llwasaki fc Tsuribel (|2008l ) found a family of 
self-similar solutions describing the condensation of radia- 
tive gas layer assuming a simple power-law cooling rate. 
llwasaki fc Tsuribel (|2009l ) investigated the linear stability 
of the self-similar solutions, and they suggested that the 
condensing layer will fragment in various scales as long as 
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the transverse scale is larger than the Field length. Physics 
of the bistable fluid has bee n inve stigated by many au- 
thors. IZel'dovich fc Pikel'nerl (Il969l . hereafter ZP69) and 
iPenston fc Brownl (Il970h investigated the structure of the 
transition front connecting the CNM and WNM under the 
plane-parallel geometry. The thickness of the transition front 
is characterized by the Field length belo w which the TI is 
stabilized by heat conduction (|Fieldlll965h . They found that 
a static solution is obtained at a so-called saturation pres- 
sure. The transition front becomes a condensation (evapo- 
ration) front if surrounding pressure is larger (less) than the 
saturation pressure. These properties of the transition front 
depe nd on its geometry . In th e spherical symmetric geom- 
etry, iGraham fc Langerl (1 19731 ) investigated isobaric flows. 
They found that spherical clouds have a minimum size below 
which they inevitably evaporate (more general de scription 
is found in iNagashima. Kovama fc InutsukJ l2005h . Linear 
stability of the transition l ayers has been investigated by 
llnoue. Inutsuka fc Kovamal (|2006l ) in the plane-parallel ge- 



ometry, and they found that the evaporation front is unsta- 
ble against corruga tion type fluctuation, wh ile the condensa- 
tion front is stable. Istone fc Zweibell l|2009t) have considered 
the effect of magnetic field on the instability. The physical 
mechanism of the corrugatio n instability is analogous to the 
Darriues-Landau instability (jLandau fc Lifshitdl 19871 ). Tak- 
ing into acc ount magnetic field perp endicular to the normal 
of the front, IStone fc Zweibell (|2010| ) investigated linear sta- 
bility of transition layers. 
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In the above-mentioned theoretical works with respect 
to the bistable fluid, the phase transition proceeds in a 
quasi-static manner through the heat conduction. However, 
recent multidimensional numeri cal simulations show more 
a dynamical turbulent structure. iKritsuk fe Normanl (|2002h 
have investigated the non-linear development of the TI in 
the three-dimensional hydrodynamical simulations with 
periodic boundary condition without any external forcing. 
As the initial condition, they set the hot ionized medium 
(T = 2 x 10 6 K) where cooling dominates heating. The 
gas is quickly decomposed into two stable phases (CNM 
and WNM) and the intermediate unstable phase. In the 
multiphase medium, the supersonic turbulence is developed 
by conversion from the thermal energy to the kinetic 
energy. They found that turbulence decays on a dynamical 
timescale. This decaying timescale is larger than that in 
supersonic isothermal turbulence of one-phase medium 
where it is smaller or compa r able t o t he flow crossing time 
Stone. Qstriker fe Gammid 1 19981 ). IKov ma fe Inutsukal 



200 6.) have done similar calculations on a much longer 



timescale. As the initial condition, they set an unstable gas 
in thermal equilibrium state with density fluctuations. They 
found a self- sustained turbulence with velocity dispersion of 
~ 0.2 — 0.4 km s _1 for a period of at least their simulation 
time (several 10 My r) in the bistable fl ui d afte r temporal 
decaying found in IKritsuk fe Normanl (|2002t ). External 
forcing such as sh ock compression can drive stronger long- 
lived turbulence. iKovama fe Inutsukal (|2002l ) investigated 
development of the TI in a shock compressed region by 
using a two-dimensional hydrodynamical simulation. They 
found that the velocity dispersion is as large as several km 
s~\ which is larger than the sound speed of the CNM. The 
supersonic turbulent motion is maintained as long as the 
shock wave continutes to propagate and provide postshock 
gas. They suggested that this supersonic translational 
motion of the CNM is observed as the supersonic turbu- 
lence in the ISM. The cloudlets are precursor of molecular 
clouds because they are as dense as molecular clouds. After 
their work, many authors have investigated the formation 
of the C NM in shocked gase s in t w o-dimensional calcu- 
lations llAudit fe H ennebelld 120051 : iHennebelle fe Auditl 
120071 : iHeitsch et alll2005l. l200rJl. and in three-dim en sional 
calculations (ICazol. Vazquez-Semadeni fe Kim 



2005; 



IVazouez-Semadeni et all 120061 : lAudit fe Hennebelld l2010h 
The influence of the magnetic fiel d on the development of 
the TI has been investigated by Inoue fe Inutsuka J2008l. 
I2009T) IHennebelle et all (|2008D . lHeitsch. Stone fe Hartmannl 
1 20091 ). 

Although physics of the bistable fluid is developing, the 
physical mechanism of driving turbulence is not fully un- 
derstood. As mentioned above, recent numerical simulations 
have found that the turbulence involving the TI is dynami- 
cal compared with the ZP69 picture that predicts very slow 
transition between CNM and WNM through the heat con- 
duction. Moreover, in the shock compressed region, because 
of its high pressure, the bistable fluid cannot exist. The ZP69 
picture cannot be applied directly to the TI in the shock 
compressed region. 

In this paper, as a first step to understand the dynam- 
ical turbulent structure, we develop the works of ZP69 into 
a more dynamical transition front. In the steady solutions, 
there are two parameters, the pressure of the CNM and the 



mass flux across the front. ZP69 considered solutions only 
on a l ine in the parameter space. lElphick. Regev fe Shavivl 
l|l992h found steady solutions with various mass flux. How- 
ever, they used a simple cubic function of the cooling rate 
that enables them to investigate analytically, and they as- 
sumed spatially constant pressure. We systematically search 
steady solutions in larger parameter space even in the high 
pressure range where the bistable fluid does not exist by 
using a realistic cooling rate. 

This paper is organized as follows: Basic equations and 
numerical methods are described in section [2] In section [3] 
we briefly review the solutions connecting the CNM and 
WNM derived by ZP69. In section [4] we present solutions 
describing dynamical transition layer that has larger mass 
flux than ZP69. The results are discussed in section [5] and 
are summarized in section [6] 



2 BASIC EQUATIONS AND NUMERICAL 
METHODS 

Basic equations for radiative gas under the plane-parallel 
geometry are the continuity equation, 



the momentum conservation, 
0, 



the energy equation, 



dE d_ 

dt dx 



dT 



-p£(p,T), 



(1) 



(2) 



(3) 



and the equation of state, 



P = TLpT (4) 

where E = pv 2 /2 + P/(7 — 1) is the total energy, and 
7 = 5/3 is the ratio of specific heats, k is the coefficient 
of heat conductivity, and £(p,T) is the net cooling rate 
per unit mass, 1Z = /cb/jtih is the gas constant, and tuh 
is the hydrogen mass. Throughout in this paper we assume 
a gas consists of atomic hydrogen. For the range of temper- 
atures considered, since the gas is almost neutr al, we adopt 
Hi = 2.5 x W 3 VT cm -1 K _1 s" 1 (|Parkerlll953h . In this pa- 
per, we adopt the following fitt ing formula of the net cooling 
rate l|Kovama fe Inutsukall2003 ). 

pC(p,T) = -B- (-T + -£-A(T)) erg cm" 3 s" 1 , (5) 
ran \ ?7iH / 

T = 2 x 10 -26 ergs -1 , (6) 

^ = >»M-^) + »-»- 2 ^ (-1 ) -< 7 > 

Figure[T]shows the thermal equilibrium state of this net cool- 
ing rate in the (n, P) plane. The gas is subject to the cool- 
ing (heating) above (below) the thermal equilibrium curve. 
The bistable fluid consisting of the WNM and CNM can 
exist in the pressure range of P m in < P < Pmax, where 
Pmin/fcB = 1596 K cm -3 and P max /fc B = 5012 K cm -3 . The 
saturate pressure where the static front is realized is as large 
as -Psat/fcB = 2823 K cm -3 in the plane-parallel geometry. 

In steady state (d/dt = 0), equations ([TJ and Q can 
be integrated with respect to x to give 
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We have two free parameters: P c and j. Given P c , the density 
and temperature can be obtained from equation Q and the 
equilibrium condition £(p c ,T c ) = 0. The momentum flux is 
determined by P c and j from equations (JSj and (|9]). 

The properties of differential equation (|14[) can be well 
understood as trajectories in the phase diagram (9,d9/dx). 
From equation (|14|l . one can see that the spatially uni- 
form state with thermally equilibrium (£ = 0) corresponds 
to a stationary point because dO/dx = d 2 0/dx 2 — 0. 
The topological property of the CNM at i = — oo cor- 
responds to saddle point llElphick. Regev fc Shavfvl Il992l ; 
iFerrara fc Shchekinovlll993r i~ Given P c and j, we numeri- 
cally integrate equation (|14|l from x = — oo along one of the 
eigenvectors obtained from the Jacobi matrix of equation 
(|14Jl around the stationary point of the CNM. 



Figure 1. Thermal equilibrium state of the net cooling rate 
C(p,T) in (P, n) plane. 



pv=j 
and 

pv 2 + P = M, 

respectively. The energy equation ([3]) becomes 



AT 



dP 



d 



p ^ dx dx dx 



(«f) +/ ,£(p,r) = o, 



(8) 



(9) 



(10) 



where j and M denotes the mass and momentum fluxes, 
respectively, that are spatially constant, and c p = jTZ/^—l) 
is the specific heat at constant pressure. From equations Q, 
(JSI , and ((9| , density, pressure, and velocity can be expressed 
by using temperature, the mass flux, and momentum flux as 
follows: 



2 i_ (m + ^M 2 - 41lj 2 T 



p (M - y/AP - ATZfT^j , 



2.7 



(11) 
(12) 

(13) 



For convenience, instead of T, we introduce new variable 
9 defined by d9 = ndT . Using 9 and equations (|lip - (|13[) . 
equation (|10|1 can be rewritten as 



P = - (m + ^AP - ATI j 2 T 



d 2 9 _ Kj 
dx 2 ~ 2k(9) 



7 + 1 



M d9 

+ ™s Tx +pc - (14) 



7- 1 ^JM 2 - A1lj 2 T{9) 



IShchekinov fc Ibanej 



(|200ll ) and 
llnoue. Inutsuka fc Kovamal (|2006l ) described a detailed 
numerical method to solve equation (|14|l . Since equation 
(|14[) is the second order differential equation, we need two 
boundary conditions. At x = — oo, we consider a uniform 
CNM with thermal equilibrium, £.(p c ,P c ) — 0, where 
physical variables in the CNM are denoted by the subscript 
of 'c'. The boundary conditions are given by 



9(x 



6> c , and 



d0_ 

dx 



0. 



(15) 



3 ZEL'DOVICH & PIKEL'NER SOLUTIONS 

First, we review the solutions connecting two thermal equi- 
librium states (CNM and WNM). We call the solutions 
derived by ZP69 the ZP solutions. The ZP solutions ex- 
ist only in the pres sure range of Pmj n < P c < -Pmax (see 
figure[l|. Given P c , Shchckinov & Ibancz (2001) determine 
j — jzp(Pc) as the eigenvalue problem so that the boundary 
condition d9/dx = is satisfied in the WNM at x = oo. 

For the case of a static solution, that is j = v = 0, 
equation (|14[1 can be rewritten as 



pCd9: 



upCdT = 



(16) 



(ZP69) where the subscript 'w' denotes the physical quanti- 
ties in the WNM. Equation (|16[1 shows the balance between 
cooling and heating inside the transition layer. The pressure 
satisfying equation (116[1 is called saturation pressure. 

In steady solutions with non-zero mass flux j ^ 0, inte- 
grating equation (|10p from x = — oo to x = oo, one obtains 



J = 



where q = 



pC(p,T)dx, 



(17) 



and we neglect the second term of equation (| 10|) since the 
flow is subsonic and hence almost isobaric. In the satura- 
tion pressure, q = since j — 0. For P > P sa t, q is positive 
because cooling dominates heating inside the front. Thus, 
j < from equation (|17|) . i.e. the front corresponds to 
the condensation. In contrast, for P < P sa t, the front de- 
scribes the evaporation. In the ZP solutions, the mass flux 
j is very small. Thus, the energy is mainly transferred by 
heat conduction. From equation (|10l) . one can see that the 
thickness of the front is ch aracterized by the Field length 
|Belgelman fc McKed[l990l ') defined by 



kT 



(18) 



The typical values of If is as large as 10 s pc in the 
CNM, and 10" 1 pc in the WNM. From equation (fT?)) . 
the flow speed across the transition layer can be evalu- 
ated by \v\ ~ lp/t c , wher e t c is the cooling timescale 
l|Ferrara fc Shchekinovlll993l '). Since this value is very small 
in the actual ISM, ZP69 suggested that the transition front 
is almost static. 
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Figure [2f a) represents properties of the ZP solution for 
Pc/ks = 3000 K cm" 3 . The corresponding mass flux is as 
large as j Z p = -3.2 x 10" 3 m H km s -1 (u w = 4.7 x 10~ 3 
km s _1 ). The first row of figure [2] shows trajectories of so- 
lutions with various boundary conditions in the phase di- 
agram (8,d8/dx). There are three stationary points, which 
are indicated by the filled circles. The CNM and WNM are 
denoted by 'C and 'W, respectively. As mentioned in sec- 
tionHJ the CNM and WNM correspond to saddle points. On 
the other hand, the stationary point at the intermediate 8 
corresponds to the unstable phase denoted by 'U' in figure 
2fa), and it is spiral point (|Elphick. Regev fe Shavfvlll992l : 
Ferrara fe Shchekinovll993h . The solutions starting from the 
CNM are denoted by the thick solid line. The ZP solution 
connects two saddle points 'C and 'W. The second row of 
figure Ufa) represents the temperature profile. The trajec- 
tory of the ZP solution in the (n, P) plane is shown in the 
third row of figure [2f a). 



4 DYNAMICAL TRANSITION LAYERS 

4.1 Transition Layers for P m in < P c < Pmax 

In the ZP solutions, the mass flux is determined so that the 
CNM is connected with the WNM. Here, we consider the 
condensation front (j < 0). What happen if the mass flux is 
larger than |jzp|? 

4-1.1 Classification of Solutions 

It is found that solutions for larger mass flux can be di- 
vided by four type solutions: C-U-o, C-U, C-c, and C-C. We 
describe the property of each solution below. 

C-U-o Solutions. For the case with j = 3jzp, the so- 
lution does not reaches the saddle point of the WNM but 
approaches the spiral point of the unstable phase (see the 
first row of figure [2p). From the second row of figure Efb), 
one can see the oscillation of the temperature profile that 
corresponds to the spiral motion in the phase diagram. The 
third row shows the trajectory of the solution in the (n, P) 
plane. The solution truncates at the heating-dominated re- 
gion before arriving at the WNM. The truncation point cor- 
responds to the peak of the temperature in the second row 
of figure HJb). Since the mass flux is still low, the pres- 
sure is almost constant. This type solu t ions already has 
been found in lElphick. Regev fe Shavrvl (|l992T ) under the 
isobaric approximation and a mo re simplified cooling rate. 
Ilnoue. Inutsuka fc Kovamal ((2006) also found them as the 
finite extent solutions, and they truncated the solution at 
the first peak of the temperature profile near the transition 
front. We call this type of solutions 'C-U-o' solutions, where 
'C-U' means solutions connecting CNM and unstable phase, 
and 'o' means the oscillation. 

C-U Solutions. For larger mass flux j — 25jzp, the 
topological properties of the stationary point of the unsta- 
ble phase changes from spiral to node (see the first row of 
figure [5Ji) El Thus, from the second row of figure [2fc) , there 

1 To be precise, before the C-U-o solution switches into the C- 
U solution, the spiral point already changes into the node point. 
In this paper, we categorize solutions with temperature peak as 




n [cm -3 ] 



Figure 3. Trajectories of solutions for j/jzp = l(the solid line), 
1000(the dashed line), 2000(the dotted line), 3000(the dot-dashed 
line). The CNM pressure is set to P c = 3000fc B K cm -3 . The filled 
circle shows the state at x = oo. 
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Figure 4. Trajectories of solutions for j/jzp = 300(the solid 
line), 500(the dashed line), 1000(the dotted line), 1500(the dot- 
dashed line). The CNM pressure is set to P c = 4000fc B K cm -3 . 
The filled circle shows the state at x = oo. 

is no oscillation in the temperature profile, and the CNM is 
monotonically connected with the unstable phase. Figure [3] 
shows solutions for much larger mass flux, j/jzp = 25, 1000, 
2000, and 3000. One can see that as the mass flux increases, 
the density (pressure) of the unstable phase at x = oo in- 
creases (decreases) along the equilibrium curve. Because of 
large mass flux, the pressure increases significantly toward 
the CNM. This typ e solutions with lower mas s flux already 
has been found in lElphick. Regev fc Shavivl (|l992h under 
isobaric approximation. We call this type of solutions 'C-U' 
solutions. 

C-c Solutions. Figure [4] shows the solutions for j/jzp = 
300, 500, 1000, and 1500. Here, we set P c /k B = 4000 K 
cm -3 . The corresponding mass flux is jzp = — 2.15xlO _2 mH 

the C-U-o solution instead of the difference of the topological 
property of the unstable phase. 
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(b)j=3j ZP (C-U-o) (c)j=25j ZP (C-U) 




Figure 2. Properties of solutions for j/jzp = 1(a), 3(b), and 25(c). The first row indicates the trajectories of solutions with various 
boundary conditions. The obtained solution satisfying equation J 1 5 1) are shown by the thick lines. The labels 'C, 'U', and 'W denote the 
CNM, unstable phase, and WNM. The second row shows the temperature profile in each mass flux. The third row represents trajectory 
of the obtained solution in the (n, P) plane. 



km s _1 . The solution for j = 300jzp belongs to the C-U solu- 
tion. For larger mass flux j = 500j'zp, the solution truncates 
at a certain point before arriving at the equilibrium state. 
This critical point corresponds to T = T cr it = M 2 / (AlZj 2 ) 
above which solutions do not exist (see equation (|lip - (|13| )). 
The physical variables at the critical point is denoted by us- 
ing the subscript of 'crit'. From equations ([8]) and ((9J, the 
condition T — T cr ;t is equivalent to 

|«crit| = \J fcrit/Pcrit- (19) 

Note that equation (|19[) is slightly smaller than the adia- 
batic sound speed y^-Pcrit/pcrit- These critical velocity and 
temperature are also seen in the steady state structure of the 
s hock front with heat condu ction without physical viscosity 
Fl (|Zel'dovich fc Raizeij[l967l ). At the critical point, one can 



The inclusion of physical viscosity may change the properties 
and cxistances of these solutions. 



see that equation (|14[) becomes infinity. Thus, in order to 
obtain a finite solution, d6/dx must vanish at the critical 
point. From equations (|13|) and (|14p . the pressure gradient 
becomes finite as follows: 

= ( 2 ^ 2 dT ) (20 ) 
V dx J crit \ JM 2 - AllpT dx j 

\ * / crit 

y J crit 

The pressure at the critical point can be derived from the 
momentum flux conservation. At the CNM (a: = -co), since 
the velocity is still subsonic because of its high density, the 
momentum flux becomes 

M = p c v 2 c + P c ~ P c . (21) 

On the other hand, from equation (|13|l . the pressure at the 
critical point is given by P cr it = M/2. Therefore, the pres- 
sure at the critical point can be expressed by P c as follows: 
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1000 2000 3000 4000 5000 6000 

-31 



P c /k B [K cm- 3 ] 



Figure 5. Classification diagram of obtained solutions in pa- 
rameter space (P c , j). The white dashed line corresponds to jzp- 
The solid line indicates the critical mass flux jcrit (Pc). The black 
area where \j\ < \ jzp\ is not focused in this paper. 



-Pcrit — -Pc/2. 

The density at the critical point is given by 

Pcrit = j 2 /-Pcrit ^ 2j 2 /P c . 



(22) 



(23) 



One can see that P cr it for j = 500jzp and lOOOjzp is roughly 
equal to P c /2 = 2000fcp, K cm -3 , and its density increases 
for larger mass flux as seen in equation (|23[) . Around the 
critical point, since dT/dx — d6/dx/n — 0, the trajectories 
in the (n, P) plane approach constant temperature lines of 
T = T C rit in figure [2 

We investigate the topological properties at the critical 
point in the phase diagram. The Jacobi matrix of equation 
(|14[) around the critical point becomes infinite. Thus, there 
is no solution passing the critical point smoothly. In order 
to obtain infinite extent solution, shock front is required. 
Solutions with shock front are described in Appendix [A] Wc 
call this type of solutions "C-c" solution, where "C-c" means 
solutions connecting the CNM and the critical point. 

C-C Solutions. For much larger mass flux case j = 
1500jzp, since the solution passes the equilibrium curve be- 
fore arriving at the critical point, the solution becomes the 
infinite extent solution connecting the CNM and CNM (see 
figure[4]). We call this type of solutions C-C solutions, where 
'C-C means solutions connecting the CNM and CNM. 

Figure [5] shows classification of solutions in the param- 
eter space (P c , j) for P m in < P c < Pmax- The white dashed 
line corresponds to j — jzp(Pc) on which the ZP solutions 
exist. The black area where |j| < \jzp\ is not focused in this 
paper. From figure [5] one can see that C-c solutions exist 
only for P c > 2P m i n . This is because if P c < 2P m i n , the 
solutions meet the equilibrium state before arriving at the 
critical point (see equation ()22[0 . 
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Figure 6. Temperature profiles for j = jzp(the solid line), 
Jcrit (the dashed line), 5jd.it (the dotted line), and 10j cr it(the dot- 
dashed line). The CNM pressure is set to P c /fcp, = 300 K cm" 3 . 



4-1-2 Structure of Transition Layers 

We investigate how the structure of transition layers de- 
pends on the mass flux. As mentioned in section [3j in the 
ZP solutions, the thickness of the front is characterized by 
the Field length because the mass flux is very small. As |j| 
increases, the effect of advection (the first term of equation 
1)101) ') becomes important compared with the heat conduc- 
tion (the third term of equation (|10[) ) in the energy transfer. 
We can define a critical mass flux jcrit where the effect of 
advection becomes comparable to the heat conduction as 
follows: 



Ik 

If.i 



T 

± m 

'F,m 



Jcrit 



jlc 



" pm Cm 7 



(24) 



where the subscript 'm' denotes a reference state in the so- 



lutions, £ Cim is the cooling length defined by 



Ic.jn — Cn 



(7 - l)p m £r, 



(25) 



and c m = \J "fPm/pm is the sound speed. As the reference 
state, we consider the state where \C\ has the maximum 
value for T > T c under the constant pressure P = P c . This 
state roughly gives minimum Field length and cooling length 
in the solution. Thus, the critical mass flux is a function of 
P c - The actual value of j C iit is shown by the solid line in 
figure [7] In most range of P c , jcrit lies in the region occupied 
by the C-U solutions. 

For Ij'I > I jcrit |, advection dominates heat conduction 
in energy transfer. In this case, the thickness of the front 
is characterized by the cooling length instead of the Field 
length. In typical ISM, If/Ic is roughly as small as 10 -2 . 
Therefore, as the mass flux increases, the thickness of the 
front increases. Figure [6] shows the temperature profiles for 
j = jzp, jcrit, 5 jcrit, and lOjcrit- The CNM pressure is set 
to Pc/fep = 3000 K cm~ 3 . One can see that the thickness of 
the layer increases with the mass flux. 
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Figure 7. The same as figure[5]but including the pressure range 
of P c > Pmax- The dashed line corresponds to T cr ; t = 10 4 K. 
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Figure 8. Trajectories of solutions for j/jalt = 1> 2, 5, and 20. 
The dashed line corresponds to P = P cr ; t /2. 



4.2 Transition Layers for P c > P max 

In previous section, we found solutions only for P m in < 
Pc < Pmax- However, in actual ISM, the CNM pressure can 
be greater than P max , for example, in shock-compressed re- 
gions. Our solutions can be directly extend for P c > P max . 

Figure is the same as figure [5] but it includes the 
pressure range of P c > Pmax- Figure shows that most of 
the regions are covered by the C-c solution for P c > P max . 
The area of C-U solutions spreads to the pressure range 
of P m ax < P c < 2P ma x. This is because solutions in this 
pressure range can connect with thermal equilibrium states 
before reaching at the critical point at a certain range of j. 
Figure [8] shows trajectories of the solutions for j/jait = 1, 
2, 5, and 20. The CNM pressure is set to 2.65 x 10 4 K cm" 3 . 
One can see that the solution terminates at the critical point 
whose pressure is equal to P c /2 and where density increases 
with the mass flux. 



5 DISCUSSION 

5.1 Implications 

iKovama fc Inutsukal Hood) have investigated multi- 
dimensional nonlinear development of the TI without 
external forcing. They found that the turbulence is self- 
sustained by the TI. The velocity dispersion for sufficiently 
large simulation box is as large as ~ 0.2 — 0.4 km s _1 . 
They suggested that this velocity is much larger than the 
prediction in ZP69. Our solutions exist in such a large 
mass flux. From their simulation, the mass flux becomes 
\j\ ~ 0.2 m H km s -1 (n/0.5 cm~ 3 )(n/0.4 km s" 1 ), where 
we assume that the accretion velocity to the CNM is 
comparable to this velocity dispersion. From figure [S] the 
corr esponding solutions he in t he C-U or C-U-o solutions. 

IKovama fc Inutsukal (|2002T ) have investigated the TI in 
the shock-compressed region. They found that the velocity 
dispersion is as large as several km s~ . The corresponding 
mass flux is \j\ ~ 20 mu km s ~ 1 (n/10 cm~ 3 )(i;/2 kms" 1 ) 
that is larger than the results in lKovama fc Inutsukal (2006). 
The pressure of the shock compressed region is larger than 
Pmax, indicting that the two stable phases cannot coexist. 
In shock compressed regions, the CNM clouds are embedded 
not by the stable WNM but by the unstable gas that is sup- 
plied continuously by shock compression. The C-c solutions 
in this pressure range (P > Pmax) may describe dynamical 
condensation from the unstable gas to the CNM in shock 
compression regions. 

5.2 Curvature Effect 

Our soluti ons are derived by assuming th e plan e-parallel 
geometry. Nagashi ma. Kovama fc Inutsukal (|2005t ) investi- 
gated the curvature effect of the transition front by assuming 
the quasi-steady approximation. They consider the cylindri- 
cal and spherical CNM cloud at the centre. In contrast to 
the plane parallel solutions, the solutions have three param- 
eters, the pressure of the CNM, the velocity of the front with 
respect to the centre, and the cloud radius. The second pa- 
rameter corresponds to the mass flux in our solutions. They 
solved equations as the eigen- and boundary-value problem 
so that the solutions connect the CNM and WNM. Thus, the 
velocity of the front f/(Pc, R c ) is a function of P c and cloud 
radius R c . The curvature effect enhance the heat flux by 
the thermal conduction from the CNM to the WNM. Thus, 
smaller cloud that has larger curvature tends to evaporate 
quickly. On the other hand, cold gas that has a negative cur- 
vature tends to gain mass by the condensation. Thus, the 
evolution of CNM cloud strongly depends on its shape. 

By using the quasi-steady approximation, we can easily 
take into account curvature effect in our solutions. It is ex- 
pected that solutions with larger evaporation rate connect 
the unstable phase and WNM while solutions with larger 
condensation rate connect the CNM and unstable phase. 



6 SUMMARY 

In this paper, we have investigated steady condensation so- 
lutions of phase transition layers in large parameter space 
(P c , j) much larger than previous works. We summarise our 
results as follows: 
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• In the pressure range where the three thermal equilib- 
rium phase can coexist under a constant pressure (P m m < 
P c < -Pmax), we find solutions that connects the CNM and 
the unstable phase. The solutions can be classified into four 
type solutions, C-U-o, C-U, C-c, and C-C solutions in order 
of the mass flux. The C-U-o and C-U solutions connect the 
CNM and the unstable phase with and without the oscilla- 
tion of physical quantities in the low-density side. The C-c 
solution is truncated at the maximum temperature above 
which there are no steady solutions. Combination of a C-c 
solution and a shock front provides an infinite extent solu- 
tion. The C-C solutions connect the CNM and CNM. 

• We have also found steady solutions in the pressure 
range (P c > P max ) where only CNM can be in equilibrium. 
This pressure range is realized in the shock compressed re- 
gion. In this parameter space (P c ,j), most of solutions be- 
long to the C-c solution. 

• We derived a critical mass flux j cr it above which the 
effect of advection becomes important compared with the 
heat conduction in energy transfer. For large mass flux, 
b'l > I Jcrit | , the thickness of the front is characterized by 
the cooling length instead of the Field length, so that the 
thickness becomes larger. 
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APPENDIX A: SOLUTIONS WITH SHOCK 
FRONT 

In this Appendix, we present steady solutions with shock 
front. Here, we assume that preshock gas is in thermal equi- 
librium state. The physical variables in preshock gas are 
denoted by using the subscript of "E". The position of the 
shock front x = x s h is determined so that the Rankinc- 
Hugoniot relation is satisfied between the solutions and the 
thermal equilibrium state by the following method. Given 
x — x s h, the specific energy flux is given by 

iPsh 



1 2 _ 7+1 2 

(j-i) Psh + 2 Vsh ~ 2( 7 -l) C *' 



(Al) 



where the subscript of "sh" denotes physical variables at 
x = x sil , and is the effective sound speed. The simple 
relation among v B h, «e, and c* is given by 



VshVE 



2 



(A2) 



From equation (|A2j) . one can get ve, and the Mach num- 
ber M is obtained. From the Rankine-Hugoniot relations, 
the physical variables in the preshock gas (pe and Te) are 
obtained. In general, the preshock gas is not in thermal equi- 
librium state. The position of the shock front x s h is deter- 
mined iteratively until £(pe,Te) = is satisfied by using 
the bisection method. 

As examples, we consider two cases of j = 2 and 24j cr i t 
for Pc/k-s = 3.114 x 10 4 K cm~ 3 . For j = 2j c case, the 
corresponding physical quantities of the preshock gas are 
n E = 0.57 cm -3 and M = 2.17. For j = 24j c , they are 
7iE = 75 cm -3 and M. = 2.23. Thus, the preshock gases 
belong to the WNM for j = 2j c and CNM for j = 24j crit . 
The upper panels of figures lA"Tk and I A lb show the density, 
temperature, and pressure distribution for j = 2j cr it and 
j — 24 j cr it, respectively. The lower panels of figure |A"T1 show 
the trajectories of obtained solutions in the (n, P) plane. 
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In addition, we perform one-dimensional hydrodynam- 
ical simulation in order to confirm the steady state solu- 
tions with shock front are realized. We consider colliding 
wall problem where the gas flows collide at x = with 
\ve — v c \ from x = ±00. Two shock wave propagates from 
x = outward. Shock heated gas quickly cools until the gas 
reaches the thermal equilibrium state. After that, the cold 
layer forms and its thickness grows by gas accretion. The 
thick gray line in each panel of figure [A"Tl indicate the results 
of one-dimensional simulations in the growing phase of the 
cold layer. One can see that the results of one-dimensional 
simulations are well described by the steady-state solutions. 
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Figure Al. Upper panels: Density (solid line), temperature (dashed line), and pressure (dotted line) profiles of steady state solutions 
at which P c = 3.113 X f0 4 K cm 3 for (a)j = 2j cr it and (b)j = 24j cr it ■ Lower panels: Trajectories of steady state solutions in the (n, P) 
plane. The thick gray line in each panel shows the results of one-dimensional hydrodynamical simulation. 
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